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The physical phenomenon of interest is high-speed gas dynamics 
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nearly-discontinuous changes in diffusion 

flowfield properties 




Aerodynamics 
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integrated heat load into a vehicle. 



Governing Equations 
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mass, and P is the pressure. 


The viscous stress tensor r and the heat flux vector q are defined as 
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Equation (4) with (7) is known as Stokes’ hypothesis for a Newtonian 
fluid [2], 


Governing Equations 
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Godunov’s theorem [4] is particularly relevant for numerical methods 
applied to high-speed gas dynamics: 
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scheme is necessarily nonlinear 
which have important implications going forth on the inteiplay 
between upwinding, shock capturing, and solution limiting. 
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Standard Galerkin procedure & SUPG-upwinding applied to (9) to produce 



Upwinding is required to stabilize convection-dominated flows, 
compressible flows t S upg is generally diagonal [3], 



SUPG stabilization does not yield monotone solutions. Additional treatment 
is needed to prevent spurious oscillations in regions of shockwaves. 
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Inviscid Flux Discretization 
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Inviscid Flux Discretization 
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contrasts to the typical approach in which 



Consider a one-dimensional, inviscid, normal shock at Mach 5. For this 
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which implies that pu, pir + P, and puH are all constant. 



Inviscid Flux Discretization 
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The semidiscrete weak form in equation (12) is discretized in time using a backwards finite 
difference scheme. Both first and second-order accurate in time schemes may be derived from 
Taylor series expansions in time about U„+i : 
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Time Discretization 
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Time Discretization 
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Implicit Solution Strategies 
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I have used matrix & matrix-free GMRES 
with full ILU-0 preconditioning - 
linearization is important 



An experimental test program was conducted in 1998 by France’s Office National d‘Etudes et 
de Recherches Aerospatiales (ONERA) to investigate shock-shock interactions produced by an 
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Static temperature contours 
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Type IV Shock Interaction 
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Type IV Shock Interaction 



0 (degrees) 



Type IV Shock Interaction 



O'dO 



Hypersonic flow over a missile nose tip with a forward facing cavity is considered. This 
configuration, shown schematically below, has been observed to exhibit transient flowfield 
response in both experimental investigations and numerical simulations. [9, 10] The flowfield 
response characteristics are largely driven by the cavity length-to-diameter ratio (L/D). 
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Forward-Facing Cavity 
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Forward-Facing Cavity 
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were conducted to assess time convergence. 
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AEDC Sharp Double Cone 



investigate the presence of 
vibrational nonequilibrium in the 
freestream [131 





AEDC Shaip Double Cone 
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Steady states, runs 2893 and 2894 
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High speed schlieren, run 2890 
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Computed schlieren, run 2890 
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Noise Characterization [14] 
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Variation of Pitot Pressure Fluctuation 
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Preliminary Results - Flowfield 




Preliminary Results - Surface Pressure 
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